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Abstract 



O 

O ■ We study the solutions to the wave equation in a two-dimensional tube of 

unit width comprised of two straight regions connected by a region of constant 
^ , curvature. We introduce a numerical method which permits high accuracy 

H ' at high curvature. We determine the bound state energies as well as the 

transmission and reflection matrices, T and TZ and focus on the nature of the 
resonances which occur in the vicinity of channel thresholds. We explore the 
dependence of these solutions on the curvature of the tube and angle of the 
bend and discuss several limiting cases where our numerical results confirm 
analytic predictions. 
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I. INTRODUCTION 



The problem of a quantum particle confined to a two-dimensional tube consisting of 
two straight regions connected by a region of arbitrary curvature but constant width, d, 
has aroused interest in recent years because of its applications in both the physics of small 
devices and electromagnetic waveguides. Related geometries such as crossed wires 

^ and elbows |^ have also been studied extensively. The primary focus has been on the 
existence of unanticipated bound states and on the "adiabatic" limit where the radius of 
curvature is always much greater than the width of the tube. 0] Goldstone and Jaffe proved 
that such a tube always has a bound state, with energy below the continuum threshold at 

In this paper, we study the specific case where the curved region is an arc of a circle. 
This case can be solved numerically with relative ease, providing a laboratory in which to 
explore phenomena we believe to be quite general but difficult to demonstrate in arbitrary 
geometries. We take the tube to have unit width. We use the arc length, s, along the outer 
boundary of the tube as one coordinate, and the distance y, along the normal from the outer 
boundary toward the center of the circle as the other coordinate. We mark the boundaries 
of the curved region as ±so. The angle subtended by the arc of the curved region is defined 
to be 6', and it satisfies the relation 9 = 2kso, where k is the curvature (here a constant). 
Figure |I| shows two examples of such a tube. Notice that we allow ourselves to consider tubes 
with 6 > 71 although such configurations require excursions from strictly planar geometry. 
Also, we choose units such that tv^ /2m = 1 making E = k'^ dimensionless. 



s 




FIG. 1. An infinite two-dimensional tube of unit width comprised of two straight regions con- 
nected by a region of constant cm'vature k. (a) 9 < ir. (b) 9 > 2tt. 

Scattering in a bent tube is a relatively simple example of a multichannel problem. 
Incoming and outgoing waves are labelled by an integer, n, the number of non-trivial nodes 
in the transverse wavefunction, smrniy. The scattering is described by transmission and 
refiection matrices Tmn and TZmn whose dimension grows as each successive channel threshold 
(at En = n^TT^) is passed. In Ref. ^ it was pointed out that the same argument which 
demonstrates the existence of a bound state, also proves that there is a "quasibound" state 
just below each channel threshold. The quasibound state below the N^^ channel threshold 
would be stable and normalizable if all lower channels (n < N) were artifically closed. Such 
quasibound states should appear as resonances in the open channel scattering amplitudes. 

Our aim in considering this solvable special case is to explore the systematics of the 



2 



bound states as they depend on k and 6, and to elucidate the character of the channel 
threshold resonances. Other workers have studied quantum wires with circular bends. Lent 
studied scattering in tubes with 9 = 7i/2 and showed that even for the largest curvature 
tractible there is very little reflection except at values of the energy just above the threshold. 
P) Sols and Macucci have confirmed this result and discovered the presence of resonances 
associated with quasi-bound states just below propagation thresholds. Moreover, they 
found that at large angles, more than one bound state may develop, and that the binding 
energy of these bound states increases as the curvature or angle is increased. 

In §11, we establish a formalism for solving for both bound and scattering states in 
tubes with circular bends. We begin traditionally - by representing the solution in the 
curved region in terms of Bessel functions and imposing the standard continuity conditions 
at s = ±so. Normally one would match the open channels to incoming and outgoing waves 
and the closed channels to falling exponentials. The matching conditions generate a linear 
algebra problem, which leads eventually to an eigenvalue problem (for the bound states) 
and a solution for 7^„ and TZmn (for the scattering problem). However, we find that this 
approach leads to serious convergence problems at high curvature. At this point we abandon 
the traditional approach and formuate a novel method (as far as we are aware) where we 
match the interior wavefunction to both rising and falling exponentials in closed channels. 
We then seek solutions by minimizing the coefficients of the rising exponentials. We argue 
that this approach should converge more rapidly and reliably. We believe this method 
has wider application than the problem at hand. Also in this section we derive a small k 
approximation to the transmission amplitude in regions far from thresholds and resonances. 

In §111 we study the convergence of our method of calculating both scattering and bound 
states. We show that it is more convergent and more efficient than the traditional method 
(of allowing only falling exponentials in closed channels). 

In §IV, we present results for the bound state eigenenergies and the scattering coefficients 
as functions of curvature and angle. We show that at fixed curvature, further bound states 
appear with increasing 6. As ^ — > oo, so does the number of bound states. The binding 
energy of the ground state approaches a limit determined by the root of a simple Bessel 
equation a.s 6 oo. For k = 1 the limit is the square of the first zero of Jq (~ (2.405)^). 
The higher bound state energies approach the same limit as the ground state, but more 
slowly. We also find the angle at which the first excited state appears, as a function of 
K, and give a theoretical explanation why this angle approaches 27r as k approaches zero. 
We also plot the scattering coefficients T and TZ as functions of the energy and verify 
the presence of resonances associated with quasi-bound states just below the propagation 
threshold for each channel. Finally, we explore the behavior of TZ and T in the complex 
plane. We explain why the elastic transmission amplitude Tu drops to zero at the energy 
of the quasibound state just below the second threshold, E = 47r^, an effect which is echoed 
in other components of T at higher thresholds. At the same time we verify the predictions 
of the small k approximation for low curvature. 



3 



II. THEORY 



In this section we summarize the basic formahsm common to both the bound and scatter- 
ing state problems. We then speciahze to the bound state case and finally, to the scattering 
case. 



A. General Formalism 

We are interested in the solutions to the wave equation, — V^^/' = k^tp within a "tube" 
coordinatized by s and y. ip vanishes on the boundaries, y = and 1, and we require 
that ip and and its normal derivative be continuous at s = ±Sq where the curvature of 
the tube jumps abruptly from zero to k. We treat the curved and straight regions at the 
same time, taking k = for the straight region. The coordinate differential is given by 
dx = sds{l — K,y) + ydy, and the differential area element is d'^x = dsdy{l — ny). Thus, the 
wave equation becomes 



d 



1 



Ky ds"^ 

We now make the standard change of variable, ip 



1 



dy"^ 
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— Ky)'^ ds^ 4(1 — Ky) 



(1 - Ky)k''ip. (1) 
■p so that eq. (|1|) becomes 

i = A;V, (2) 



If we define r = 1 — /t?/, and separate variables we recognize the radial equation as Bessel's 
equation, with solutions of the form. 



ipiis.r) = —Si{s)yi{r) 

1 Y 



where 



yi{r) = Vr 



Y, 



kr 



J. 



\ ^ J 



k) 



and 



Si{s) 



cos KUiS 

sin KVis 



(3) 



(4) 



(5) 



Here Ji, and Yi, are Bessel and Neumann functions respectively. The linear combination of 
Bessel functions in eq. satisfies 3^ = 0at?/ = l(r = 1 — k). The order, vi, must be chosen 
to satisfy 3^ = at y = (r = 1). vi may be either real or imaginary, as we discuss below. 
The label / = 1, 2, . . . denotes the number of non-trivial nodes (/ — 1) in y{r). 

To study the character of vi, we substitute eqs. (^) and (|^) back into the wave equation. 
We find that the momentum satisfies the relation, 



(6) 
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where is a positive term that grows with /. Thus, for fixed k"^ only a finite number of 
eigenvalues, uf, can be positive. For any k"^ greater than vr^ there is at least one positive u"^ - 
a fact which is closely related to the existence of a bound state. As k increases, more positive 
z/^ values appear - associated with the channel threshold resonances discussed below. Let 
I < i correspond to real u and I > i give imaginary z/, where we define (3 = iv. Then, 



y/r \ sm Ki/is J y/r \ smh/t/J/S J ^ ^ 



We have added the subscript "<" to ip to denote the solution inside the circular region. 

Depending on k"^, the exterior solutions, \s\ > sq, either oscillate or vary exponentially 
with s. The general exterior solution is of the form, 

oo 

^> = ^(fo^e"^™ + c^e^'"^'-'")) sinrnTry. (8) 

m=l 

where 



7« 



Vm%2 _ p. (9) 



When k'^ > m^vr^, 7^ becomes imaginary and the solutions become oscillatory. In that case 
we define 7™, = ikm with km > 0. 

We now apply continuity conditions at s = sq, recalling that in the interior region the 
normal derivative of is 

s-VV^< = -|-V^<, (10) 
r OS 



and we find 



bm + Cm = 2y aiSi{so)M, 



ml 

1=1 



-bm + Cm = —yai—^\^^^^Nmi (11) 

Im ^_-|^ Wii 

where M^i = dysinmiry^ and N^i = dy sin rmry^^. Eqs. (|TTp give an infinite 
dimensional matrix linear algebra problem to solve for the energies of bound states or the 
iS-matrix for scattering states. In practice we must truncate the sum over internal and 
external channels. We superpose L internal solutions and attempt to match to M external 
channels. 

When k'^ < N'^n'^ the {cm} with m > N are coefficients of exponentially rising solutions. 
Because these solutions are not normalizable, standard practice is to dismiss them as un- 
physical and set {cm} for m > N to zero from the beginning. This leads to an eigenvalue 
condition for k which has unique solutions only when L = M, i.e. one is forced to use 
channel spaces of the same dimension both inside and outside the region of curvature. As 
the dimension of this channel space (L = M) is increased, the accuracy of the solution 
should improve. After all, it corresponds to taking more terms in the approximation to 
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the wavefunction. However, we have found that this straightforward approach fails at high 
curvature (k > 0.6). As the channel space dimension is increased, the solution does not 
seem to converge. We believe that this is due to poor matching at the boundary s = Sq: We 
expect that an interior function with n nodes should match to a superposition of exterior 
functions dominated by terms with n, rail, and n ±2, etc. nodes. However, the L = M 
restriction forbids the introduction of any exterior functions with higher frequencies than 
the interior functions, so that the n + 1, n + 2 node terms are not present. Clearly it would 
be advantageous to allow M > L. However, a simple count of the number of degrees of 
freedom shows that the there is no solution under these circumstances. 

To solve this problem, we have implemented a novel approach: we retain the rising 
exponentials and seek first the set of {ai} and next the value of k for which Ylm=i'^m 
minimized. As we shall see, this approach removes the L = M restriction, allowing us to 
take as many terms in either sum as we like. Of course, in the limit L, M — > oo, {cm} 
in the closed channels for all physically interesting solutions, whether bound or scattering. 



B. Bound state formalism 

Bound states correspond to solutions of Eqs. ([TT]) with A;^ < vr^. Our strategy will be to 
allow the {cm} to be non-zero and seek solutions which minimize J2m=i ^m- We anticipate 
that the parity of the ground state will be even, so we choose internal solutions which are 
superpositions of cos kus or cosh KjSs. If there are several bound states, this symmetric 
ansatz will find the even ones. A corresponding ansatz odd under s ^ —s finds the odd 
states. Even and odd bound states alternate as a function of the energy. 

We use the method of Lagrange multipliers to minimize ^ subject to the constraint 
that Xl^i '^f — 1- There is some arbitrariness in this prescription. The {a;} must be 
constrained in some way, otherwise the trivial solution {a;} = will follow. The sensitivity 
of our method to the form of the constraint is discussed in the following section. We define 

M L 

A({a,}, k)^Y. ^-({«'» - ^£ - 1) (12) 

m=l 1=1 

and seek to minimize A first with respect to the {ai} and then with respect to k. The 
conditions |^ = 0, lead to an algebraic problem, with a solution {a^{k)}. Substituting back 
into eq. (|T^, A reduces to a function of k, whose minima we locate numerically. 
First we combine eqs. (Ill]) to obtain function of the {a/}, 

L 

Cm = '^Gmiai, where (13) 
1=1 

Gml = M^iSiiso) + —KmiS'iiso). (14) 

7m 



Next we substitute for Cm in eq. (|T2D, and differentiate with respect to ai to obtain. 
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''^^Eijttj = \ai, where 



M 

'^Ij = QmlGmj (15) 

m=l 

This is an L dimensional, matrix eigenvalue problem. The external channel space di- 
mension, M, appears in the definition of S, but does not complicate eq. (|T5|) . Since S is 
a positive definite matrix, it has real, positive eigenvalues, A^(A;), \^{k) etc., ordered from 
smallest to largest, and associated eigenvectors, a|, of, etc. Each eigenvector gives a sta- 
tionary point of the functional A. Normalizing Ylii=ii.^\Y = 1 as required by the constraint 
it is easy to see that 

M 

A({aKfc)},fc) = J]c^({ai}) = V(fc), (16) 

m=l 

SO the smallest eigenvalue gives the minimum A(fc) at fixed k. To find the best approximation 
to the bound states for fixed L and M, we search for a value of < vr that minimizes the 
lowest eigenvalue, A^(fc). The magnitude of \^{k) measures the squared sum of coefficients 
of rising exponentials, and is a measure of the accuracy of our solution. There is always at 
least one minimum of \^{k) for k < n — a consequence of the general theorem of Ref. [0. 
For small values of Sq, there is only one minimum at A; = A;i, corresponding to a ground state 
energy of kf. As sq is increased (at fixed k), there will be multiple minima of X^{k) which we 
denote by ki, k^, k^, etc., with kp < n. {ki corresponds to the ground state, k^ to the second 
excited state, and so on.) The odd-numbered excited states {p even) are antisymmetric in s 
and are obtained from an antisymmetric ansatz. 



We must set the number of terms in the sums in eqs. (p^Sf ) to values such that the 
computations are tractable while still retaining reasonable accuracy. Studies of convergence 
and accuracy are presented in §111. We emphasize, though, that the convergence is not a 
trivial issue for this problem: the default assumption, L = M, fails to converge as L — »• cx) 
for interesting values of k, forcing us to take a different approach. Anticipating the results 
of §111, we note that the results presented in §IV were obtained with L = 4 and M = 10, 
which generated accuracy to about five decimal places. 



C. Formalism for scattering states 

We now extend our formalism to handle scattering solutions. At an energy E = k"^, the 
"longitudinal" momentum in the n^^ channel is kn = Vk'^ — n^vr^. If A^^vr^ < k"^ < (A^-|-l)^7r^, 
then A^ channels are "open" : for n < N, kn is real and wavefunctions in the external regions 
are oscillatory. For n > N, kn is imaginary. These channels are "closed": we define ikn = jn] 
the solutions are exponential as in the bound state problem. Since internal and external 
channel sums cover different ranges we will use different indices as follows: i, j, k, and I 
label internal (curved region) solutions; and m, n, p, q, and r label external (straight region) 
channels. A^ will always denote the number of open channels. 
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An incident wave in the g*'* channel will couple to all channels. However, the solutions in 
the closed channels die exponentially, so that at a distance (s ±C)o), we see reflection and 
transmission in only the open channels. We wish to solve for these TZ and T coefficients and 
explore their behavior as functions of the channel and energy of the incoming wave as well 
as the parameters of the tube. Unlike the bound state problem, this is not an eigenvalue 
problem. There are many solutions (in fact 2N) for each value of the energy. Our aim is to 
find the solutions and parametrize them in terms of the matrices TZ and T. 

Since our internal solutions are naturally constructed as either symmetric or antisym- 
metric (in s), it is convenient to do likewise with the external solutions. We consider the 
symmetric case explicitly — from which we will obtain T -\- TZ. The antisymmetric case 
proceeds along identical lines and yields T — TZ. The solutions are labelled by the indices q, 
the channel of the incoming wave, and p, the channel under consideration. The symmetric 
solution is the sum of the left- and right- incoming solutions: 

s ^ ±00 : al{s) = -i=(<5p^e-^^-l^l + (T + 7^)p,e+*'=''l^l) (17) 



Unitarity requires {T + 7Zy {T + 7Z) = 1, and from the antisymmetric solution, {T — TZy {T — 
TZ) = 1, or 

r+r = 7^"^7^ = i, 
r^n + n^T = o. 

(18) 

These constraints can be summarized by the requirement that the matrix 

7^ T 
T TZ 



S=[ll] (19) 



be unitary, S'^S = 1. This is a special case (for Hamiltonians symmetric under s — >■ —s) of 
the general definition of the iS-matrix, 

where TZ and T are reflection and transmission coefficients for waves incident from the left 
and TZ' and T' are the analogous coefficients for right incident waves. Time reversal invari- 
ance requires = S, which requires = T and TZ^ = TZ for symmetric Hamiltonians. 

m 

To solve for the matrices T and TZ, we return to the basis wavefunctions for the bound 
state problem, modified to accommodate the open channels and the factor of y/ik. Note 
that this factor is present only in the open channels where n < N. Because our solution is 
symmetric, we take Si{s) to be cos kuis, 
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^ I 



^ 1 



ra=l 
A/ 



(21) 



n=Af+l 



By comparing Eqs. (0) and ( plD we can identify T + 7^ in terms of 6 and c for open 
channels, 



d = (r + 7^)p,e 



ikpso 



(22) 



for p,q < N. By analysis similar to the bound state problem, Q continuity at s = sq shows 
that 



and 



where 



1=1 



1=1 



(23) 



(24) 



/Ik:Siiso)M^i + -^S'iiso)N^i n<N 
^ 5Kso)M„z + :^5Kso)iV„; n>N 



(25) 



and 



tknSi{so)Mni - -^S'i{so)N^i n<N 



n> N 



(26) 



It is convenient to redefine the internal wavefunction Si in the case where f/ is imaginary. Instead 
of cosh kPi we use ^[e~'^/^'^^o-s) _|_ g-K/3!(so+s)j_ 'pj^jg amounts to multiplying by a constant e"'^^'*", 
which will not change our result but keeps the numbers tractably small when sq is large. We 
must now examine where vi is actually imaginary. For real v the s-dependence of takes the form 
cosKi^s = ^(e*"^'^* + e~^'^^'^). By strict analogy with the straight regions, we would then expect 
real v for open channels / < N and imaginary u for / > A^. However, the energy thresholds for 
propagation are actually lower in the curved region |^, so that we will find, for example, 1^2 i^eal 
beginning not at k"^ = Att"^, but at k"^ = An'^ — €2, where €2 is a number that we can find numerically 
because we expect = at the threshold. 
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As in the bound state case, tradition would dictate setting the coefficients of growing 
exponentials in closed channels to zero ab initio. We seek to improve convergence and 
accuracy by instead keeping the growing exponentials and minimizing their contributions. 
We have found that the best convergence for high k is obtained by minimizing the ratio of 
the growing contributions relative to the falling contributions, i.e. minimizing Yln>N l*^"!^ 
subject to the constraint that J2n>N ~ ^■ 

Because scattering boundary conditions are complex {viz. eq. (P^D), we must be more 
careful about the reality properties of the variables we consider. Certainly c„ and 6„ are 
complex in open channels. Likewise, the matrices Ti. and Q are complex, although they are 
clearly real for n > N. Our manipulations will be simplified considerably by using matrix 
notation. Even though Q and Ti are not square matrices, there should be no ambiguity in 
the definition of quantities such as Q^TC, etc.. We define a matrix projection onto the closed 
channels: Vmn = 0, except Vmm = 1 when M > m > N . Note that P is an M dimensional 
matrix in the external channel space. 

For an incoming wave in the g*^ channel (suppressing the q superscript), we can once again 
define the projection onto rising exponentials, and implement the constraint Yln>N l^nP ~ ^ 
by a lagrange multiplier, 

A(a, k) = a^Ea - \{a}m - 1) (27) 

where S = Q^VQ and = 1-OVH. Note that H and Vt are both real, symmetric matrices in 
the internal (L-dimensional) space. In practical calculations, L is small (compared to the 
external channel space dimension, M), so these matrices are relatively easy to manipulate. 

It is convenient to work in the basis of eigenvectors of VL. Let {v'^,k = 1,L} be the 
orthonormal eigenvectors of Q corresponding to (real, positive) eigenvalues {uJk}, 

v''^-if = 6ki. (28) 
We rewrite a in the basis of the {■u^}, using y/uJ^ as a metric, 

^ 1 

S=y2^Vkv'. (29) 

Inserting this definition back into eq. (^Tf ), we obtain 

A{{f]},k) =f]^Ef]- \{f]^f]-l). (30) 



where 



'-'^^ - ^—_v'^Ev' (31) 



Now, the condition ^ = becomes 



^a^ 

dA 



O^Ef] = \f], (32) 
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an L dimensional eigenvalue problem with solutions {ff'andA'^, forg = As in the 

bound state case, the eigenvalue A'' is the value of A at the stationary point. Small A'^ there- 
fore correspond to solutions to the wave equation with negligible projection on exponentially 
growing solutions in the closed channels and are therefore physically interesting. 

The number of physically significant solutions to this eigenvalue problem depends on 
the energy, fc^. Below the lowest threshold, fc^ < vr^, we do not expect (and did not find) 
any solutions with small A except at the energies of bound states. When one channel is 
open (tt^ < k"^ < 47r^), there should be a single symmetric solution to the wave equation for 
each energy. So we expect to find one eigenvalue, A^, much smaller than the rest. When 
a second channel opens, there are now two symmetric solutions, so we expect to find two 
small eigenvalues, A^ and A^. We conclude that we must expand the vector a as a linear 
combination of the A^ eigenvectors associated with the lowest A^ eigenvalues when A^ channels 
are open. We now expand the vector a in terms of the lowest A^ eigenvectors ff of S: 

N L ^ 

a = Y.drY,^vlv\ (33) 

r=l k=l ^'^^ 

where dr is the expansion coefficient and the sum over k corresponds to the expansion in 
terms of t]^ as defined in eq. (^). Finally, we define the L ^ N matrix 

■^^r = (34) 

so that we may write compactly 

d = Md. (35) 

We are now prepared to return to our identifications of 6^ and from eq. ([2^). Inserting 
our above expression for a, into eq. (0), we find that 

d'n = Qnqe''^'^ (36) 
where Q is an A^ x A^ matrix defined by 

Q = KM. (37) 

Substituting eqs. ([35| ) and (^) back into eq. (p3D and recalling our identification with (T+TZ) 
in eq. (0), we obtain at last 

(T + n)p, = {GTWQ-^ (38) 

Note that in the end the matrix manipulations which are calculationally intensive — finding 
eigenvalues, eigenvectors and inverses — are only carried out on matrices of dimension L. 
The external channel index, M, can therefore be taken very much greater than L without 
increasing the computation time significantly. 

The antisymmetric scattering states can be analyzed in the same fashion, beginning with 
the parameterization of the scattering states, 

s ±oo : alis) = -^{6„e-'''-\^\ + (T - 7^)p,e+^'=-l^l) (39) 

obtained by subtracting the left- and right-incoming wave solutions. This allows us to solve 
for the quantity (T — TZ)pq, which we then combine with eq. ( |38D to yield the desired matrices 
T and 71. 



11 



D. Small curvature approximation for the transmission amplitude 



When the curvature, k, is small, the solutions to our problem simplify somewhat. There 
is a quite well developed theory of the adiabatic approximation for a quantum wire with 
small and slowly varying curvature — k <^ 1 and ^ ^ 1 — but it does not apply to a 
circular region joining two straight regions because n jumps abruptly to zero at sq. Also, 
most of the results are restricted to sqrt{k'^ ~^^) ~ 0. Nevertheless we can use the fact that 
K is small to simplify the calculation of the elastic scattering amplitudes Tmm and TZmm in 
regions away from resonances and bound states. The result is 

7;„„^[lH-0(.= )]exp(^,^ 

n^mm = y sin(2fc^so - k^9/2) [1 + 0{k)] exp +~0{f kO) ) ' 

where = \/W^—t^ and these expressions are only valid when k'^^ k,"^. 

The derivation of these results begins from eq. (0), which we expand in a Fourier series, 
(P{s,y) = En/n(^) sinnTT?/, 

Dmnifn + ^fn) + kifm = 0, (41) 

where 

^ f\ sin miry sin mry 

Since -Dmn is of order k, when m 7^ n, it is easy to see that the effect of the n^^ channel on 
a wave incident in the m*^ channel begins at order k"^. This confirms our decision to work 
only to order k. To this order Dmm = 1 + which makes this system unusual, since the 
perturbation appears in lowest order only in the kinetic term. In the end this will lead to 
the result that an arbitrarily small curvature can result in a large scattering phase provided 
the product usq (x 6 is large. 

Working to order k we can ignore the "potential" with the result that fm oc 

exp -^{iqms), where = A;^(l — We now construct symmetric and antisymmetric 
combinations of these interior solutions, and match them and their normal derivatives to 
the appropriate external scattering states eqs. ([17|) and (pQ)). Since channel coupling is 



of higher order in k, Tmm ± Tlmm are both unitary, so we parametrize them by phases, 
Tmm ± 'Tlmm = exp(2i(5™). The result of the matching calculation is 

tan(A;„so + 5™ ) = (1 ± ^ + 0{k^)) tang^so. (43) 

Note that qmSf) and kmS^ differ by —kmHSo/2 which need not be small even when n 
because 2kso = 6 can be held fixed in this limit. So the leading effect on Tmm and TZmm arises 
from this accumulating phase along the bending section. The prefactor (1 ± |) determines 
the modulus of Tmm and Rmm,- When we solve for \Tmm\, we find that effects of order k 
cancel, which they must do since \TZmm\ ~ ^ and jT^mP + I'^mmP must give unity to order 
K^. This completes our derivation. The results, eqs. (^OD, give a useful approximation to the 
elastic amplitudes over much of their range when n <^1. 
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III. COMPUTATION AND CONVERGENCE 



The straightforward approach to solving the wave equation by matching interior to ex- 
terior solutions failed for physically interesting values of curvature and angle, forcing us to 
develop a different method of calculation. We are not aware of any studies of this method in 
the literature, and therefore present some details of our studies of convergence here. Readers 
not interested in calculational methods should skip directly to the results presented in §IV. 

A. Convergence of calculations for the bound state energies 

The traditional method of matching L independent internal solutions to L independent 
external solutions does not converge as L increases at large values of k. This is illustrated 
graphically in Figure ^ where we plot the calculated ground state energy ki as a function 
of L for K = 0.8 and 6 = it. We have attributed this failure to converge to the absence of 
small wavelength (in y) states on the outside, necessary to accurately expand the internal 
states. In order to demonstrate that this is indeed the case, we have studied the convergence 
as a function of L and M. Table | shows the effect on the calculated ground state energy 
ki of varying the channel spaces L and M in the case where k = 0.8 and 6 = n. The best 
accuracy for the least amount of computation time is obtained by setting M > L. We can 
see from the lower left corner of the table that setting M < L does not produce favorable 
results, because as higher-frequency interior functions enter the picture, we require even 
more high-frequency exterior functions for a good match at the boundary s = sq. 

3.1 16| ' ' ' 1 



3.114- 



3.112- 




■ 5 10 15 20 

FIG. 2. The eigenvalue ki plotted as a function of the channel space in the case L = M, for 
K = 0.8,6 = IT. As the channel space is increased, the solution starts to converge but then begins 
to oscillate wildly. 

We notice a few general trends in the convergence of ki. As expected, increasing M at 
fixed L improves the accuracy of the result. Note that the larger the value of L, the larger 
M must be taken to capture this improvement. Note also that once the improvement has 
occured, there is little to be gained by further increasing M. So, for example, L = 5 improves 
dramatically as M increases from 5 to 10, and improves little thereafter. In contrast L = 9 
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L 


M=5 


M=10 


M=30 




M=50 


1 


2 9837e-03 


2 9834e-03 


2 9834e-03 




2 Q834e-03 


2 


5 5981e-05 


5 4785e-05 


5 4635e-05 

W^T.WWWVx' WW 




5 4635e-05 

Wi^EWWWV^ WW 


3 


8 3369e-05 


7 2461e-05 


7 2181e-05 




7 2180e-05 

1 t^-Lwwv,^ WW 


4 


6 0521e-06 


9 3112e-06 

• w -±- -i- ^ w w 


9 1729e-06 




9 1717e-06 

w ■ -L 1 J.- 1 w w 


5 


7.4863e-05 


1.0638e-05 


9.9300e-06 




9.9287e-06 


6 


2.6707e-03 


1.6007e-06 


2.2744e-06 




2.2728e-06 


7 


7.6036e-03 


4.1820e-06 


2.3406e-06 




2.3390e-06 


8 


8.3919e-04 


3.6752e-06 


5.7210e-07 




5.7000e-07 


9 


4.7389e-03 


9.8571e-06 


5.7590e-07 




5.7370e-07 


10 


9.7321e-04 


6.9994e-05 


2.9000e-09 




O.OOOOe+00 




TABLE I. Tabulation 


of \ki{L,M) 


- A^i(10,50)| for K = 0.8, 


9 = IT. 





improves through M = 30, and httle thereafter. Also, for fixed M, the convergence worsens 
as L is increased past the value of M, as we can see clearly in the M = 5 case. Increasing 
the number of interior basis functions should improve the accuracy of the representation 
of the interior wavefunction, but apparently this is useless when the corresponding higher- 
frequency terms needed to match to the outside wavefunction are omitted, as is the case 
when M < L. As can be seen from the table, our parameters of L = 4 and M = 10 yield 
accuracy to about one part in 10^. 

We also note that as the number of interior basis functions L is increased, the value of 
^(^i) = decreases, again provided that M > L. Table |I| illustrates this convergence 

for the same parameters k = 0.8 and 6 = tt. This is to be expected; as the bound state 
approximation is improved, the contributions of the exponentially rising terms should ap- 
proach zero. Although the value of A(fci) does not decrease for fixed L and increasing M, 
we suspect that this is a result of the increased number of terms in the sum Ylm which 
offsets the expected decrease in rising exponential contributions. 

For the case L = M, A{ki) deviates from the trend in that it is many orders of magnitude 
smaller than for the other cases. This can be explained by recalling the traditional approach 
discussed in §11. A, where the coefficients of the rising exponential exterior functions are set 
to zero from the beginning. The resulting set of equations has solutions only when L = M. 
With our approach, we thus expect to find a solution where the sum of these coefficients 
A{k) is zero (limited only by machine precision) when L = M. 

B. Convergence of different methods for scattering states 

In our study of scattering we explored further the choice of constraint in the Lagrange 
multiplier problem. One obvious choice is simply to constrain our solution to match to the 
asymptotic conditions in a channel of interest, namely that the equations for 6^ (eq. p^)) 
must be satisfied in the open channels. Our equation for A would then require a Lagrange 
multiplier for each separate channel, leading to 

M N 

A= ^ {G*^ia*i){Grnjaj) -^\n{Hniai-bn). (44) 

m=N+l n=l 
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I, 


M=5 


M=10 


M=30 


M=50 


1 

X 




9 94sq7fip n"^ 


9 94R89^p (T^ 


9 94RS98p (T^ 


9 






1 90379Qp-n4 


1 90381 4p-n4 


■? 






4 9fi1 791 p-f)^ 


4 9fii Qn8p-n^ 


4 


7 4499fiQp-nfi 


8 Qn74nip-nfi 


8 Q^n7fiip-nfi 

O. C/tJW 1 \J ±\j \J\J 




5 


2.754173e-19 


4.409149e-06 


4.686405e-06 


4.688083e-06 


6 


2.209942e-17 


1.565973e-06 


1.632015e-06 


1.6332936-06 


7 


1.440247e-09 


7.985669e-07 


1.019960e-06 


1.021530e-06 


8 


1.770510e-05 


3.232110e-07 


4.602362e-07 


4.613563e-07 


9 


1.573734e-04 


1.585803e-07 


3.173988e-07 


3.188958e-07 


10 


7.881601e-02 


2.996393e-19 


1.676304e-07 


1.685153e-07 




TABLE II. A(/c"i) 


for K = 0.8, 6 = TT with varying 


I channel space 


L and M. 



We note that the incoming channel label, g, introduced in §11. C, which should appear on 
A, ttj, and has been suppressed for clarity. 

Extremizing this expression for A yields a set of equations from which we can determine 
ai and A„ in a straightforward manner. However, we have found that this method has poor 
convergence for large k. Next we considered the same condition as used for the bound states, 
namely eq. (|T^ with the extension that the coefficients ai may now be complex. While this 
second method converges for significantly higher values of n than the previous method, we 
get even slightly better convergence with the procedure outlined in §11. C, in which we use 
the constraint ^ = 1. 

Table illustrates the convergence of |7^ii| for different L and M in the case k = 6.28, 
K = 0.9, and 6 = tt /2 for the constraints Hniai = bn and ^ = 1. This is a region — near 
the second threshold at 2n and at large k — where convergence is problematic. By examining 
|7?.ii(L, M) —7^11 (14, 24)1 for the two choices of constraint, we see that convergence is better 
for the second method. This is primarily due to the fact that for a given value of L, increasing 
M continues to improve accuracy for a longer range in the second method. For example, 
in the case L = 10 ,M = 12, the two methods are comparable, but when we reach L = 10, 
M = 22, the quantity |7^ii(L,M) — 7^(14, 24) | is a full order of magnitude smaller for the 
second method. 

As with the bound state problem, we find that we must take M > L to achieve conver- 
gence. To illustrate this further. Figure ^ plots \TZu\ for the same parameters, using the 
method outlined in §11. C, as a function of L for M = 50 and for M = L. From this, we 
can see clearly that the M = L case is unsatisfactory, as expected. Although the speed of 
convergence varies with the region of interest — it is significantly slower when the energy is 
near a resonance — we have taken parameters which give accuracy to at least a part in 100. 
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M 




A 


u 


o 


1 n 


1 9 


14 


2 


2.4008e-01 














4 


3.1798e-01 


3.7502e-01 












6 


3.2236e-01 


1.4818e-01 


1.2767e-01 










8 


o r\ o /i o r\'\ 

3.2343e-01 


1.6741e-01 


4.7922e-02 


2.4390e-01 








10 


3.2381e-01 


1.7536e-01 


4.8338e-02 


3.2088e-02 


1.3033e-01 






12 


3.2397e-01 


1.7893e-01 


5.2158e-02 


1.7015e-02 


4.8322e-02 


3.8592e-01 




14 


3.2405e-01 


1.8067e-01 


5.4878e-02 


1.5969e-02 


1.1385e-02 


1.0571e-01 


1.8553e-01 


16 


3.2410e-01 


1.8159e-01 


5.6600e-02 


1.6559e-02 


6.2273e-03 


1.7310e-02 


2.2753C-01 


18 


3.2412e-01 


1.8209e-01 


5.7671e-02 


1.7235e-02 


5.4160e-03 


4.4498e-03 


4.1930e-02 


20 


3.2414e-01 


1.8240e-01 


5.8343e-02 


1.7773e-02 


5.4195e-03 


1.9425e-03 


7.7122e-03 


22 


3.2415e-01 


1.8258e-01 


5.8771e-02 


1.8168e-02 


5.5804e-03 


1.3692e-03 


1.5033e-03 


24 


3.2415e-01 


1.8270e-01 


5.9049e-02 


1.8451e-02 


5.7499e-03 


1.2636e-03 


O.OOOOe+00 




(a) |7^ll 


(L,M) -7^( 14, 24)1 using 


the constraint Hnioj — 


6^ for all n 


< N. 


M 




4 


6 


8 


10 


1 9 


14 

J. 4: 


2 


2.4129e-01 














4 


2.3975e-01 


1.2888e-01 












6 


2.3969e-01 


4.3621e-02 


1.3154C-01 










8 


2.3968e-01 


3.1081e-02 


2.3833e-02 


1.8675e-01 








10 


2.3968e-01 


2.7993e-02 


9.4324e-03 


2.7598e-02 


2.7725e-01 






12 


2.3968e-01 


2.7046e-02 


5.8291e-03 


7.0832e-03 


4.6534e-02 


3.6196e-01 




14 


2.3968C-01 


2.6712C-02 


4.6140C-03 


2.5538e-03 


9.8018C-03 


9.1442C-02 


4.0612e-01 


16 


2.3968e-01 


2.6580e-02 


4.1339e-03 


1.1407e-03 


2.4954e-03 


1.8614C-02 


1.7307e-01 


18 


2.3968e-01 


2.6523e-02 


3.9258e-03 


5.9470e-04 


4.5970e-04 


4.3728e-03 


4.0580e-02 


20 


2.3968e-01 


2.6496e-02 


3.8295e-03 


3.5300e-04 


2.5870e-04 


7.7920e-04 


9.5638e-03 


22 


2.3968e-01 


2.6482e-02 


3.7825e-03 


2.3630e-04 


5.5810e-04 


3.6960e-04 


2.1651e-03 


24 


2.3968e-01 


2.6475e-02 


3.7584e-03 


1.7670e-04 


6.9850e-04 


8.0860e-04 


O.OOOOe+00 



(b) |7^ll(L, M) - 7^(14, 24)1 using the constraint J2 l^nP = actual method. 



TABLE III. Tabulation of |7^ll(L,M)| - |7^ll(14, 24)| for k = 6.28, k = 0.9, and 9 = tt/2 for 
two different choices of constraint. 
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2 4 6 8 10 12 14 2 4 6 8 10 12 14 

L (M=50) L (M=L) 

(a) (b) 

FIG. 3. |7^ll| vs. L for (a) M = 50 and (b) M = L. Note that, as with the bound state 
problem, setting M = L produces unsatisfactory results. 



IV. RESULTS AND DISCUSSION 

A. Eigenenergies of the bound states 

We first examine the dependence of the bound state energies on the parameter sq, the 
length of the curved region of the tube. Figure ^ plots the momenta kj corresponding to 
bound state eigenenergies against the angle of the circular arc 6 (= 2k,sq) for curvature 
K = 0.8. Eigenenergies for both symmetric states and antisymmetric states are shown. 




10 20 30 40 50 60 

e 



FIG. 4. Plot of kj vs. 9, for k = 0.8. 

The curved region introduces a term that behaves like an attractive finite square well 
potential with depth Increasing the length of this curved region both lowers the energy 

of the ground state and allows for more excited states. 
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3.1 

e=i6.5ji 

FIG. 5. Plots of A(k) vs. k for k = 0.8 and different 9. Only symmetric states are considered. 
The minima correspond to kj for odd values of j. 

Figure |^ shows plots of the lowest eigenvalue of the matrix S as defined in §11. B against 
k, for curvature k = 0.8 and different values of the angle 6. The minima on these graphs are 
candidate (symmetric) bound states provided k < tc, and correspond to points in Figure § 
along a line of fixed angle. As the angle increases, the number of eigenstates also increases. 

In the limit 6 — > oo, all the kj approach the same value, which we will denote by koo- 
The value of fcoo is that which corresponds to u = for any given curvature. To see this, 
note that in the limit ^ — oo, the external straight region becomes negligible (for bound 
states), and the internal curved region becomes very long. In this limit it becomes possible to 
accomodate wavefunctions with any number of nodes (in s) without introducing significant 
energy. The longitudinal wavefunction, cos kus, (the antisymmetric case can be handled 
analogously) can oscillate any number of times over the range —sq < s < so contributing 
only K^z/^ to the energy. By choosing u very small this contribution to the energy can be 
made negligible. A minute tuning of u allows one to match the wavefunction to a falling 
exponential at s = Sq. Thus as Sq ^ C)0 there are an infinite number of bound states 
accumulating at the value of k corresponding to z/ = 0. At any fixed very large Sq there are 
also many other bound states at energies between k^o and vr. The value of k^o is determined 
by solving the transcendental equation. 



Yo { M Jo (^(1 



0. 



(45) 



This provides us with koo{K.). As k ^ 1, Yq becomes singular in eq. (^51). A careful analysis 
of the behavior at small 1 — k shows that the equation reduces to Jo(^oo('^ = 1)) = 0? which 
gives fcoo(l) = 2.405.... fTTJ 

Figure D plots the solution of eq. (|45|) , koo{n) against ln(l — k). As k — (ln(l — k) — >• 0), 
the tube is straight and so u = corresponds to k = n, as expected. And as k — 1, 
koo — ^ 2.405 ... as expected. It should be noted that the amount of binding increases 
dramatically fact which necessitated our choice of ln(l — k) as the independent 

variable in Fig. 10. 



18 



2.4L 



-35 



-30 



-25 



-20 -1 5 
ln(1-K) 



-10 



FIG. 6. The value of k^, function of ln(l — k). 

Next we turn our attention to determining the angle at which the first (antisymmetric) 
bound excited state appears. We label the angles at which new bound states appear as 9i, 
02, etc. where 9i corresponds to the ground state and 9j (j > 1) to the (j — 1)*'' excited 
state. For any k, 6i is zero, since there exists a bound state for any tube with a curved 
region, no matter how short. 



6.9 
6.8 
6.7 
6.6 
6.5 
6.4 
6.3 
6.2. 



0.2 



0.4 0.6 

K 



FIG. 7. The angle 9 at which the first antisymmetric bound state appears as a function of 
curvature k. As k ^ 0, 6 ^ 2tt. 

Figure shows a plot of 62{i^)- As can be seen from the graph, 62 approaches 2tt as 
K approaches zero. We can see that this must be the case by considering the one-channel 
approximation to the interior wavefunction, namely 



ijj = sinvrysinz/iKS. 



(46) 



It is sufficient to consider this approximation because in the limit that k —>■ 0, the tube is 
very close to straight and thus the contributions from higher channels are negligible. In this 
limit, the problem reduces to that of a square well potential. The Schrodinger equation then 
requires 
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2 2 ; 2 2 

Z/, K = k — TT . 

4 



(47) 



At the angle where the bound state first appears, k = tc so eq. (^) imphes z/i = i. For a 
bound state to exist, the interior wavefunction (sinz/i/ts) must be able to match to a falling 
exponential at the boundary s = sq. The smallest value of sq at which this can occur is that 
which satisfies the condition uiksq = tt/2. Recalling that Sq = 9 /2k, we then have 

We note that as n is increased, the value of 62 also increases. This seems somewhat 
counter-intuitive, as we expect from the Schrodinger equation that higher curvature implies 
greater binding, so that the length of the curved region necessary to bind an excited state 
should decrease with large k. This is indeed the case. However, we must recall that 9 
depends on k as well as sq. For a fixed value of 6, sq is much larger for small k than for 
large k. Thus, the attractive potential acts over a much longer distance for small k, and this 
effect overshadows the difference in binding energy due to differences in n alone. 



B. Scattering results 

Previous workers have noted many features of the scattering problem, and we do not have 
any qualitatively new phenomena to add to their treatments. The primary features of 
interest to us are the rapid rise of |7^i| from threshold, the steady increase of the phase 
of %nm. in regions of energy away from thresholds, and finally, the dramatic fiuctuations in 
T and TZ at energies just below thresholds. Our interest lies in understanding better the 
physical origins of these effects. These features are due to the presence of a bound state 
just below threshold, a steadily growing phase over a wide range of k, and weakly coupled 
quasibound states just below each channel threshold, respectively. The last of these provides 
an elegant example of a resonance manifesting itself as total destructive interference — a 
phenomenon familiar to particle physicists in the form of the sudden drop in the tttt nn 



S'-wave cross section just below KK threshold, due to the presence of the resonance. 112 
We begin by examining |7ii| and |7^ii| as functions of k. Figure||plots I^Lip and |7^iip, 
the measures of transmitted and reflected flux, for 1 < < 3 at flxed k = 0.9 and 9 = 7t/2. 
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0.4 



0.2 ' 

1 



1 1.5 2 2.5 3 

k/ji 

FIG. 8. ITiip (solid) and |7^llp (dashed) vs. fc/vr. Here k = 0.9, 9 = tt/2. 

Note that in the region where N = 1 (ji <k < 27r), jT^ip + |7?.iip = 1. This is required 
by unitarity and serves as a check on our numerical calculations. We have checked that our 
5-matrix remains unitary in the N x N well. As k crosses the second threshold at 

27r there is significant inelasticity into the newly open channel evidenced by the fact that 

|rnp + |7^np<l. 

We also observe that at energies just below A^^vr^, the transmission amplitude in the 
highest open channel (7/v-iAf-i) drops sharply almost to zero. We shall show that this and 
other rapid variations in T and TZ are manifestations of resonances related to the presence 
of quasibound states just below each new channel threshold. They are the energies at which 
we would find a bound state if all open channels were artificially closed. 

It is difficult to figure out what is happening by considering |T| and \TZ\ alone. We will 
learn much more by looking at the Argand diagrams for T and TZ. The general characteristics 
will be illustrated by considering plots of Tu and TZn in the complex plane for it < k < 27r 
as shown in figure 
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(a) (b) 
FIG. 9. Argand diagrams for (a) Tu and (b) TZu for the parameters k = 0.9, 9 = tt/2, 
TT < k < 2Tr. Arrows indicate direction of increasing k; circles mark equal intervals in k. 

TZu and Tu are not independent when only one channel is open. Unitarity [S'^S = 1) 
requires both + |7^iip = 1 and ReT{\7lu = 0. We can parameterize both by the 

modulus (M) and phase (0) of Tu, 



Tu = Me'^ 
TZu = iVl - M'^e'^ 



(49) 



So it suffices to discuss the behavior of Ti- According to Fig. @, Ti vanishes at k = n. As 
k increases, Ti rapidly executes a clockwise circle of radius 1/2 in the complex plane until 
\Tu\ ~ 1 at a A; only slightly larger than vr. Then Ti moves clockwise on the unit circle until 
suddenly, at a k just below 27i, Tu suddenly executes an almost complete counterclockwise 
circle of radius 1/2. Tu briefly resumes its steady phase growth, until the second channel 
opens at /c = 27r and Tu moves off the unit circle. 

First consider the behavior near k = tt. We know on general grounds that Tu = and 
7^11 = — 1 at threshold. Furthermore, we know that there is a symmetric bound state at an 
energy k'^ = kf = ir'^ — e. This must appear as a pole in T + 7?. at this energy. Since S must 
be unitary for any real /c^, Tu + TZu must be of the form. 



T 



11 



ki + ie 
T^ii = 1 + Ja{k) 



IE 



(50) 



where k\ = a/P 



and fa{k) is a smoothly varying function of k near k = tt. There is 



no nearby antisymmetric bound state, so Tu — TZu is a smooth function of k, 



Tu-n 
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near k = tt. Comparing with the threshold behavior, Tu = and TZu 
/^(tt) = and fai'^') = 1. Combining all these results we obtain. 



(51) 

— 1, we see that 
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TZn = ie^'^ sin (f) where 

e 

tan = — 
h 



(52) 



At threshold, = 7r/2, and as increases, (j) drops rapidly to zero over an interval scaled 
by e. Eqs. ( |52D parameterize a semicircle of radius 1/2 centered at Tu = 1/2, which Tu 
executes as ki grows from zero to a value greater than e. The speed at which Tu moves 
measures e, by the formula 

e = ik, - l] (53) 



For ki ^ e the conditions of the small k approximation discussed in §11 are satisfied, 
provided n is small. Thus we expect to find TZu ~ and Tu = 7^° xexp{—iki9/2). We verify 
this by comparing the rate of change of the argument of Tu with the small k prediction, 
daigTii/dki = —6/2. So in the range where the approximation is valid, we expect that the 
graph of the phase angle of T will be linear with respect to ki, with a slope equal to half 
the angle subtended by the curved region of the tube. Computing Tu for 4.0 < ki < 5.0 
with K = 0.001 and 6 = tt/2, a range and curvature for which the small k approximation 
is valid, we have found daigTu/dki ~ — 7r/4 = —9/2, which confirms that this relation is 
indeed satisfied. 

Finally we turn to the dramatic behavior of T and TZ near channel thresholds. It is clear 
from Fig. (P) that Tu resonates {i.e. it rapidly executes a counterclockwise circle in the 
complex plane [0|) at the energy where \Tu\ drops precipitously to zero. This behavior is 



indicative of a pole in the complex A;^ plane just below the real axis. The pole appears in 
Ti + T^ii because the quasibound state is symmetric in s. Tu — IZu is smooth over this 
region in k. If we denote the pole location as qi — i'ji, where qi ^ 2tt and 71 is small and 
positive (by causality). Then in the neighborhood of the pole, 

7i.+^„^ t'"'":''" +---. («) 

ki- qi + Z7i 

where . . . denotes other, smooth terms in Tu + T^ii- Combining this with a smooth form 
for Tu — IZu we see that Tu executes a counterclockwise circle of radius 1/2 in the complex 
A;-plane as k passes over the pole. Just before the vicinity of the pole, Tu lay nearly on the 
unitarity circle, \Tu\ ~ 1. The only way that it could rapidly execute a counterclockwise 
circular path consistent with unitarity is if the resonant circle is tangent to the unitarity 
circle. [Otherwise Tu would have to pass outside the unit circle at some point on its circular 
path.] Since the radius of the resonance circle is 1/2, the (tangent) resonant circle must pass 
through the origin. Thus \Tu\ must vanish in the near vicinity of the quasibound state as 
observed by Ref. 0, and shown in Fig. (^. 

The basic structure of Tu and TZu can be seen in the scattering coefficients for higher 



channels. As we observe in figure |T0| for N = 2, T12 and follow the same pattern. For 
^2 = y/k"^ — 47r2 just above zero, T executes rapid clockwise motion as a result of the pole 
located near k = n. This is followed by smooth phase growth until k reaches the quasi-bound 
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state just under the = 3 threshold, at which time it executes the rapid counterclockwise 
circle in the complex plane characteristic of an elastic resonance. 




_ii ■ 1 _ii . 1 

-1 1-1 1 

Re(T22) Re(R22) 

FIG. 10. Argand diagrams for T12 = T21, T22, Tli2 = '^21, and 7^22 for the range 2tt < k < Svr. 
The same basic structure is observed as for the coefficients in the = 1 case. 

Our conclusion from this exercise is that the behavior observed in Fig. ^is characteristic 
of scattering in bent tubes. The bound states, quasibound states and regions of smooth 
phase growth are robust properties of these systems. They do not depend specially on 
our choice of a circular geometry which facilitated our numerical calculations. Perhaps the 
most interesting observation is the appearance of violent fluctuations in transmission and 
reflection properties of bent tubes in the vicinity of channel thresholds. The fluctuations 
occur over very narrow intervals in energy (the resonances are narrow), so low resolution 
experiments would quite likely fail to detect them. On the other hand, careful experiments 
using, for example, microwave radiation in waveguides, should see these rapid fluctuations 
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in response associated with the opening of channel thresholds in which quasibound states 
occur. 
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